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The Coulomb problem for Schrodinger equation is examined, in spaces of constant curva- 
ture, Lobachevsky iJs and Riemann 6*3 models, on the base of generalized parabolic coordi- 
. nates. In contrast to the hyperbolic case, in spherical space S3 such parabolic coordinates 

^jQ' turn to be complex- valued, with additional constraint on them. The technique of the use 

' of such real and complex coordinates in two space models within the method of separation 

of variables in Schrodinger equation with Kepler potential is developed in detail; the energy 
CO . spectra and corresponding wave functions for bound states have been constructed in explicit 

form for both spaces; connections with Runge-Lenz operators in both curved space models 
Qh' ^-re described. 

^ . 

, I. INTRODUCTION 

i Quantum mechanics had been started with the theory of the hydrogen atom, so considering 

' the Quantum mechanics in Riemannian spaces it is first natural step to turn to just this system. 

, A common quantum-mechanical hydrogen atom model is based materially on the assumption of 

00 

' the Euclidean character of the physical 3-space geometry. In this context, natural questions arise: 

what in such a model description is determined by this assumption, and which changes will be 
entailed by allowing for other spatial geometries: for instance. Lobachevsky's H3, Riemann's 
53, or de Sitter geometry. The question is of fundamental significance, even beyond its possible 
experimental testing. 

Firstly, the hydrogen atom in 3-dimensional space of constant positive curvature 5*3 was con- 
sidered by Schrodinger He had been studied the so-called factorization method in quantum 
mechanics; in particular, application of this techniques to discrete part of the energy spectrum 
for hydrogen atom had been elaborated. An idea was to modify the basic atom system so as to 
cover all the energy spectrum including the region E > as well. However, mere placing of the 
atom system inside a finite box in order to make the whole energy spectrum discrete did not seem 
attractive, so Schrodinger had placed the atom into the curved background of the Riemann space 
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model S3. Due its compactness, the spherical Riemann model may simulate the effect of the finite 
box - see Schrodinger |]| and Stevenson 
In spherical coordinates of 5*3 



dx^ +sm^x{dO'^ + sin^ 



the Schrodinger Hamiltonian in dimensionless units has the form 



H 



2 ^ dx"" dx"' tan x ' 

p is a curvature radius, a unite for length; M is a mass of the electron; fi? /Mp^ is taken as a unit 
for energy; e = stands for a Coulomb interaction constant; the sign at e/tan x corresponds 

to the attracting Coulomb force. The energy spectrum is entire discrete and given by 



■2;? + 2(" 



n 



1,2,3,... 



Hydrogen atom in the Lobachevsky space was considered firstly by Infeld and Shild 



a 



H 



dr = dx + cosh^xidO + sin'' d(j) 



2 tanh x 

Energy spectrum contains a discrete and continuous parts. The number of discrete levels is finite, 
they are specified by 



— < e < (- - e 
2 - - ^2 ' 



p2 1 

7T (n 

2n2 2 ^ 



n 



1,2, 3,..., TV. 



In the region e > — e) the energy spectrum is continuous. 

Thus, the models of the hydrogen atom in Euclid, Riemann, and Lobachevsky spaces sig- 
nificantly differ from each other, which is the result of differences in three spatial geometries: 
-E'3, Hs, S3. To present time, we see a plenty of investigations on this matter: 
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II. COORDINATES IN SPACE H3 AND ^3 

In Euchdean 3-dimension space £"3 there exist 11 coordinate systems 
complete separation of variables in the Helmholtz equation 



481], allowing for the 



g°'^{x) stands for the metric tensor of space E3 specified for curvilinear coordinates (x^,x^,x^). 
Solution of the same problem for spaces of constant positive and negative curvature, Riemannian 
53 and Lobachevsky H3 models was given by Olevsky in [49]. It was established that there exist 
34 such coordinate systems for hyperbolic space H^, whereas in the case of spherical model S3 the 
number of those systems is only 6. Extension of the analysis to the complex sphere was given by 
Kalnins and Miller in 



5l| ]. Also see 



53-55 



551, 



56|. 



The above asymmetry between H3 and ^3 may be seen as even more strange if one calls the 
known relations of these models through the analytical continuation 



H3 
S3 
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Xq — Xi 



x\ 



^3 



+ 1 , 



(1) 



the curvature radius R is taken as a unit for length. 

The asymmetry of the models H3 and S3 with respect to coordinate systems finds its logical 
corollary when turning to the study of the quantum mechanical model for a hydrogen atom on the 
background of a curved space. In particular, an additional degeneracy like in the case of flat space 
was observed, which presumes existence of a hidden symmetry in the (curved space) problem. In 
J-lD], the symmetry operators accounting for such additional degeneracy in Kepler problem on 
curved space background ware found for both model H3 and 5*3, and an analog of the conventional 
Runge - Lenz vector in the flat space was constructed. 

Connection between the Runge - Lenz operator A in the quantum Kepler problem and parabolic 
coordinates in Euclidean space is well known: by solving the Schrodinger equation in these coor- 
dinates the eigenfunctions of the third component ^3 arise [13]. Analogous situation exists in the 
hyperbolic space H3 but not iri in the spherical 5*3. In the Lobachewsky space, among 34 coordi- 
nates established by Olevsky j49| one may select one special case of parabolic system of coordinates 
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in H3, in which the Schrodinger equation ahows the separation of variables and the wave functions 
arisen turn out to be eigenfunctions of the operator B = + L^. Among six coordinate systems 
mentioned in [4] an analog of parabolic coordinates is not encountered. 

If one looks at 34 and 6 systems in H3 and 6*3 respectively, one can note that all six ones from 
5*3 have their counterparts in H3. The main purpose of the present paper consists is the search 
of some counterparts of remaining 34 — 6 = 28 systems. It turns out that such 28 systems in S3 
can be constructed, but they should be complex- valued; to preserve real nature of the geometrical 
space one must impose additional restrictions including complex conjugation. 

In particular, the complex analog for parabolic coordinates in space of the positive curvature 
53 can be introduced and used in studying the quantum mechanical Kepler problem in this space 



(this possibility was partly studied before in [14 1 
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4^) 



Let us start with the following fact: from the the metrics in Lobachevsky space (note x ^ 
[0,+oo)) 



(2) 



by means of the change X ~^ > sinh x ~^ ^ sin x one can obtain the corresponding metrics of the 
Riemannian space (note that x G [0)'''"]) 



- [ dx^ + sin^ X (f^^*^ + sin^ 



(3) 



This simple observation on ^^3 — 5*3 connection leads us to interesting consequences. Indeed, let 
us compare, for instance, wave functions and spectra for hydrogen atom in spaces of negative and 
constant curvature 



5(x) = sinh'xexp 



n-l-l 



n 



X 



xF ( - + / + 1, / - n + 1, 2/ + 2; 1 - e"^^ 

„2 



1 



2n^ 



{n' - 1) ; 



(4) 



5(x) = sinxexp ( i(n - / - 1) - - ) x 



-i- + / + 1, / - n + 1, 2/ + 2; 1 - e 
n 



1/2 
2n2 2 ^ 



(5) 
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quantity [Mh?' / B?) provides us with natural unit for energy, e is a dimensionless parameter char- 
acterizing intensity of the Coulomb interaction. One may readily note that these two solutions 
turn into each other at the following formal replacement 

X — ^«X,e — y -ie , e — > - e . (6) 

This example indicates that the relation between and 6*3 reflected by substitution x ~^ 
is meaningful. In the context of the described above situation with coordinate systems in and 
5*3, let us make use of this correspondence {x — > i x) as follows. 

Let in the Lobachevsky space ^^3 be chosen a coordinate system {pi, p2, ps) (one of those 
34 found by Olevsky), then as a first step one has to establish connection of such a system with 
spherical one: 

Pk = fkix, 9, <t>) , (7) 

and a second step is to introduce a corresponding coordinate system in the space Ss through the 
formal change x — ^ ^ X- 

Pk = fkiix, 0, (j)) . (8) 

With help of this prescription one can determine 34 coordinate systems in space 6*3 in comparison 
to six ones given in 4^. It turns out that 28 new (added) coordinate systems are complex- valued 
and therefore additional restrictions should be imposed which involve complex conjugation. All 
these extra coordinate systems permit the full separation of variables in the Helmholtz equation 
on the sphere 53. 

Below, only one example of such coordinates, analog of the parabolic ones in space H-^, will be 
examined in detail and applied to the study of the quantum-mechanical Kepler problem on the 
sphere 5*3. 

III. PARABOLIC COORDINATES IN SPACE MODELS ^3, iJg 

In Q, 

the following coordinate system (the case XXV) in Lobachevsky space had been given 
= A{p,-a){p,-bY + A{p,-a){p,-bY '^P^ - - - "^^^3 ' 
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where {pi, p2, ps) are connected with the four quasi-Cartesian coordinates (xq, xi, X2, X3) 

- {Xof + (Xif + ix2f + (X3)2 = -1 , Xo > +1 



by the formulas 



With the notation 



and a = +1, 6 = 0, eq. (|T0|) gives 



X2 

— = tan [{a — b)p3] , b < pi < a < p2, 

X\ 

xj + xl^xj-x^^ix^-x^^^ 
Pi - a Pi-b [pi- by 



xl + xl = a'^ , X3 - xo = U , X3 + X0 = V 



(T^ + UV = — 1 , xi = a cos p3 , X2 = cr sin p^ , 

UV UV U^ 
+ + ^ = 0, - + + — = 



Pl-1 ' Pl ' pj ' P2-I P2 P2 

Combining two last equations results in 

pi - I P2-I pi P2 

{-^-^)'y' + {pi-P2) uv = o. 
pi - I P2-1 

From this, taking into account cr^ + UV = —1, we arrive at 

U P1P2 



V P1P2 - pi- P2 

and further 



, UV = P1P2 - pi- P2 



(10) 



jr2 TT P^P'^ - Pi- P2 

U = P1P2 , V = u . 

P1P2 

Thus, for U, V, a we have found (the Lobachevsky space is realized on the branch xq > +1, so that 
(X3 - xo) < ) 

TT I ^r I P\^ P2- P1P2 

(7 = X3 - xo = - , y = X3 + xo = -= , 

^JP\P2 

a = V-1 -UV = V-(l - Pi)(l - P2) ■ 
Explicit formulas relating pi,p2,pj, with Cartesian coordinates (xo,x;) look as 



xi = \/-(l - pi){l - P2) cos P3 



X2 = V-(l - pi){l - P2) sinp3 , (11) 

Pl+ P2- 2/91P2 Pi + P2 

X3 = , xq = ; 
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and the inverse formulas are 

Xo -X3 Xq- X3 

Pi = : , P2 = , 

xo + x xo-x ^^2) 

P3 = arctan— , x = + Xn + Xo. 

Xl \/ 1 ^ o 

Now, instead of the introduced pi,P2,P3 one can define other coordinates which behave simply 
in the limit ii — t- oo (the curvature vanishes). Such a limiting procedure for spherical coordinates 
of the hyperbolic space with metrics 

dl^ = [ dx^ + sinh^ X {(le^ + sin^ 9 dcj)^) ] 

going over into spherical ones of the flat space 

lim (px) = f 1 liin {p sinh x) = r . (13) 

Eliminating xq as follows 

qi = — = , , qi = tanh y n; , ni = (sin cos 6, sin 9 sin 6, cos 9) 

we note that when R ^ oo the coordinates qi reduce to 

lim (pqi) = lim {p tanh x ^l) = r ni . (14) 

p— >oo p— >oo 

So, to have coordinates with known and understandable behavior in the limit i? — )■ oo we define 
new coordinates ti,t2,(j) 

ti = 1 - pi = — — , t2 = l- P2 = -. , 9 = P3 = arctan — ; (15) 

1 + q 1-q qi 

in the limit of the flat space they provide us with the known parabolic coordinates (^, ??, (p) 

lim (pti) = z + r = ^ , lim (pt2) = z — r = —rj . (16) 

p— >oo p— >-oo 

The metrics ([9]) in coordinates {ti,t2,4>) takes the form 

,,2 *i ~ *2 ,,2 , h-h 2 . . 1,2 
'^^ = At,{l-t,f '^'' + 4t2(l-t2)^ "'^'^^"^ ' 

0<ti<l, t2<0, 0<(?i<27r. (17) 

Now, with the help of the rules ([7]) and ([8]), one can define corresponding parabolic coordinates 
ti,t2 on the sphere S^. To this end, coordinates (ti,t2) in ^3 must be expressed in terms of 
spherical ones (x, 0) 

tanh Y , — tanh Y 

ti= 1 + COS0 I , t2= 1-COS0 — — — ^ 18 

1 + tanh X 1 — tanh x 
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from whence we get defining relations for corresponding coordinates in ^3 

i tan Y , ^-^ —i tan y , ^ 

Take special notice that (ti and in ([19]) are complex-valued expressed through two real (X)^)- 
The inverse formulas are readily found 

1 + cos6' = + — ) , 1 - cos6l = - — ) , 
iq iq 

ti+t2 — 2tit2 . ti — 12 , . 

cose = , iq = - — . (20) 

11—12 I — ti — 12 

So defined parametrization of ^3 by coordinates ti,t2 can be additionally detailed by the for- 
mulas 



ti = (l + cose)(^(x), t2 = (l-cos0)^*(x) , 

vr 

ip{x) = sinxexp - x)] ■ 
Prom (12ip one can derive the relationship between ti and t2 

ti ti h{i-t2) 



(21) 



which are equivalent to 



^ , tit2 = (tlt2)* , (22) 



t*i - -h ' *2 - ' '•^^^ 



or 



:i-<.) = (i-<.)(-|), (i-*o = (i-y(-|). 



(24) 



Its existence may evidently be referred to the real nature of the space S^. In particular, one 
consequence is: if — ^ 1 — then t2 — ^ 1 + 0, and so on. 

In the following, so defined coordinates (ti, ^2) 4') are called parabolic coordinates on the sphere 
53. In the limit of the flat space, they reduce to the ordinary parabolic coordinates 

lim {—ipti) = ^ , lim {—ipt2) = —77 . (25) 

p— >oo p— >oo 

Now, we transform the metrics of the space 6*3 in spherical coordinates 

df = -dx^ - sin^ X (de^ + sin^ 6 dc/)"^) , 
to complex parabolic ti, t2, 4'- As a first step, with the help of 



sin 9 = tit2 2 — ' ^^'^ 



1 + ^2 q2 

^— , sm x= 
q^ 1 + 
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we obtain 

sin^ X sin^ 9 dcp^ = h t2 d(t)^ . (26) 

As a second step, we have 

('^^ = cose) = — -—^ [d{ ) , 

sm^ 6 11*2 1 + q h — 12 

and further 

sin^ xidOf = 7-^-^ ^ , , ,4 [ t2{t2 - 1) dh - h{h - 1) dt2 f , 



or 



sin2 X {dOf = 4(i_^^p(i_^^)2^^^^ [ t2 it2 - 1) dh - *i (*i - 1) dt2 f • (27) 



Taking into account relations 



and 

2 1 (2-ti-*2)^ 

COS X- - 



l + tan2x 4(l-*i)(l-*2) ' 
we get 

f..2 [{l-t2)dh-{l-h)dt2? 

^^^^ " 4(l-*i)2(l-*2)^ • ^'^^ 

Therefore, for the metrics in parabohc coordinates in S3 we have arrived at the form 

Formally, this formula differs from its counterpart in the space H3 only by presence of (—1) in the 
expression for dP. 

The more clarity in the complex coordinates may be achieved if one determines them in terms 
of Cartesian coordinates 

i yi = v'-*i*2 cos ^ , i J/2 = \/-*i*2 sin 4> , 
. _ *i + *2 - 2*1*2 _ 2 - *i - *2 _ (30) 

2V(l-*i)(l-*2) ' ^°"2V(l-*i)(l-*2) ' 

they evidently obey an identity (below the notation y = + + 2/3 is used) 



y? + yi + yi + 2/0 = (y + ^yo)(y - iyo) = i 
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Inverse to (1301) formulas are 



ti = {y + ys) {y + iyo) = « e*" , 
t2 = {y- yz){y - iyo) = h e"*" , 

^1*2 = ?/i + ^2 1 tan (/) = — . 

yi 

Two complex coordinates are detailed by 

a = y + y3, b = y -ys , ab = yf + yj , 
cos a = y , sin a = yo ; 

from this it follows 



cos a = y 



a + b a — b 



ys ■ 



2 ' 2 

The domain for variables (a, b) can ne illustrated by the Fig. 1. 
b 




(31) 



(32) 



Fig. 1. The domain for {a,b). 

Note that associated points in the sphere {yo,yk) and {—yo, —yk) are parameterized according to 

(+yo,+yfc) =^ ih, t2, (p) , 

(-yo,-2/fc) =^ (i2,ti,0 + vr) , 



(33) 



or 



(34) 



i+yo,+yk) =^ {a,b,(j)) , 
{-yo,-yk) =^ {b,a,(t) + -iT) . 
Let us describe peculiarities of parametrization of S3 by the variables (a, 6, a, </>). The first that 
is a closed line 



2/1 = 0, ?/2 = , 2/3 + 2/0 = 1; 



(35) 
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in this case, ab = 0, and to this hne there corresponds only a part of the boundary in Fig. 1, a = 
and 6 = (and the coordinate <p is "dumb") 



71" TT 

y3>0 , ti = 2^3(2/3 + iyo) , t2 = , ae[ - 2 ' + 2 ^ ' 
y3 < , ti = , t2 = -2y3(j/3 - iyo) , a G [ +| ] . 



(36) 



Now, let us consider another closed line 

y3 = 0, yo = 0, yl+yl = l; (37) 

characterized by 

a = b= ^Jyl + yl = 1 , e±^" = l±i0,a + 6 = 2. (38) 

To this line there correspond a single point (1, 1) on the boundary (see Fig. 1), and the coordinate 
^ is not "dump" now 

yi = cos (p , 2/2 = sin 4> . (39) 

There exists one other peculiar region to which there corresponds the line a + 6 = 2 on the 
boundary (see Fig. 1). Indeed, let 

yo = , yl + yl + yl = l, (40) 

this sphere is parameterized in accordance with 

sin a = , cos a = 1 ^> a + b = 2 , 4> E [0, 2-it] , 

^ ^ (41) 

yi = +\/ a{2 - a) cos (p , 2/2 = +v a{2 - a) smcp , j/3 = a - 1. 

Here one can introduce the variable a — 1 = cos 9, then 

yi = + sin 9 cos (p , 2/2 = + sin 9 sin (p , y^ = cos 9 . 

IV. SEPARATION OF VARIABLES IN PARABOLIC COORDINATES, THE MODELS 

H3 AND 53 

Now let us turn to a Coulomb problem. General expression for Schrodinger Hamiltonian 
1 1 d ^ d e ^ . ti — t2 

gives 
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S3 

1 — It d , ,9 1 — to 9 / ^ d 

h-t2dti^ ^'dti t2-tidt2 ^ 'dt2 

^ _ • 2 - tl - t2 

2tit2 9cA2 ti-t2 ■ ^ ""^ 

Transition to the Hamiltonian in the model is achieved by performing two formal replacements: 
e =^ ie, and H =^ —H. Thus we obtain 

l-h d , , d l-t2 d , , d 

H = -2 — ti(l - tl 2 — 12 1 - t2 7— 

tl-i2 5ti ^ '^9ti t2-tiat2'^ '^9t2 

^J_^ _ ^ 2-t,-t2 

2tit2 9(/>2 tl -t2 ■ ^ ^ 

Now, let us separate the variables in the Schrodinger equation in coordinates {ti,t2,4>) - first 
let us specify the case of the model S3 

^(ti,t2,0) = /i(ti)/2(t2)e^'"^ (45) 

from i7 ^' = it follows 



2(1 -tl) d d 2(1 -t2) d d 

17^ ir,''^^- wj'^^' -i^-tT dr2 ''^^ - dr2 

hh - ^e l^ll^ f,f, = e M2 . (46) 
ztii;2 tl — t2 

and then multiplying ()i6|) by (ti — t2)/2 /1/2 , one derives 

1 /- ^ d ,^ ^ d ^ m? ie e 

-^(l-'^'^*^<'-'^)^A+=^^+ft. + |t. + (t,-'=.) = (47) 

where two separation constants ki and k2 obey an additional constraint 

k^-k2 = -ie. (48) 

Thus, the problem in 5*3 consists in solving the system 
S3 



(l-t2);f t2(l-t2);f /2 + ( 

at2 dt2 



2 

-ie — € 



t2 



m 
4t7 



m 



+ h)fi = o, 



(49) 



^ + ^2)/2=0. 
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In the model H^, we obtain (separation constants obey the identity ki — k2 = e) 

N d , . d „ , e + e m , , „ 

Solutions of (|49p and ()50p are searched in the form 

/l = (1 - ti)"^ Si(ti) , /2 = t^^ (1 - ts)"^ S2{t2) . (51) 

It suffice to consider the case (I49p : transition from 6*3 to is realized trough the changes 

53 ^ , e ^ -e , +ie . 

The first equation in (j49]) gives (to obtain analogous result for /2, is suffices to change the index 1 
by 2 and the parameter e by — e) 

hil - h) S'l + S[ [ 2a(l - fi)26ti + (1 - 2ti) ] 
aiiai - 1) ( ^ - 1 ) - 2ai6i + 61(61 - 1) ( 1 - 1 ) + ai ( 1 - 2 ) - 61 (2 ^ 



ie — e ( 1 \ / I 1 \ , 1 



+— ^ . - 1 - ^ - + + k 



=1- 



5i(ti) = 0. (52) 



Terms proportional to tj^^ and (1 — ti)^^ can be eliminated by imposing additional restriction: 



ie — e m? 



which results in 



af- — = 0, bi + — _ + /ci = 0; (53) 



ti(l - ti)S'-l + S;[(2ai + 1) - (2ai + 26i + 2)ti] 

ie — e 



ai(ai + 1) + 2ai6i + 61(61 + 1) + 



Si = 0. 



2 

This means that S{ti) is expressed in terms of hyper geometric functions Si{ti) = F(ai, /3i, 71; ti) 
with parameters obeying to 

71 = 2ai + 1 , ai + /3i + 1 = 2ai + 26i + 2 , 

ZC — 6 

ai/3i = ai(ai + 1) + 2ai6i + 61(61 + 1) + , 

from whence it follows 



1 /I e-ie 

Qi = ai + 61 + - + y - + 



1 I e — ie 

/3i = ai + 6i + -- + , 71 = 2ai + l . (54) 
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Let us summarize: 
for 53 

Si = F{ai , /3i, 71; ti) , S2 = F{a2, (^2, 72; ^2) 

2 2 i2 e — ie m? o e + «e m? 

«i = — . ^2 = — ; 61 = + — -ki, 62 = + — - A;2 



1 I e-ie , 1 /le + ie 

Qi = ai + 61 + - + y - + , 02 = 02 + 02 + - + Y - + 



1 I e-ie „ , 1 /I e + ie 

/3i = ai + 6i + --y- + ^— , ^2 = 02 + 62 + --^- + ^— , 

71 = 2ai + 1 , 72 = 2a2 + 1 , ki - k2 = -ie; 



and for 



/l=tr /2 = t2Ml-t2)'^52 

= F(qi, /3i, 71; ti) , S2 = F(q2, (^2, 72; ^2) 

2 7 2 — e + e o — e — e 

X' "2 = — , ^1 = ^— + — -^!' ^2 = ^— + — -^2 



1 /I -e + e , 1 I -e-e 

ai = ai + 61 + - + Y - + — - — , 02 = 02 + 62 + 2 + Y 4 + — 2 — 



(55) 



(56) 



1 1 -e + e ^ ,1/1 -e-e 

/3i = ai + 61 + - - ^ - + , /32 = a2 + 62 + - - Y 4 + ' 

71 = 2ai + 1 , 72 = 2a2 + 1 , ki - k2 = e . 
V. THE HYDROGEN ATOM IN H3, BOUND STATES 

The task consists in separating from all solutions found in Section IV those which describe 
possible bound states. To treat this problem we will need some additional details in parametrization 
of the model H3 by coordinates ti,t2,(j) (see ([IT]) ) 



X3 



Xl = V — tlt2 COS I 
t2 — 2tit2 



X0 + X l + q 
-00 <t2<0 



t2) 

< ti < 1 



Xq 



t2 



X2 = V—'tit2 sincj) 

2-ti-t2 
v/(l-tl)((l-t2) 

X3-X qs-q 



XQ 



X 



1 



tan (h = — = — 



X2 
Xl 



q2 

qi 



e [0, 2tt] . 



(57) 



(58) 
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Besides, we will use relations 



Qi = ;r cos( 

^ 2-ti-t2 



V^hh ^/{^^tl)0^^^ . , ti+t2-2tit2 
12 = sm (p , q3 = . 

2 — ^1—^2 2 — ti — t2 

Note that the origin is characterized by 

ti=0, t2 = <^ qi = , 92 = , 93 = ; 
and the whole boundary for the domain G{ti,t2) is detailed as follows (see Fig. 2) 

(a) ti = 0, t2 G (-00,0] =^ 

t2 



qi = , q2 = , q3 = -q - , 

A — 12 

t2 -oo , qs = -q ^ -1 ; 
(b) , t2 G (-00,0] =^ 

1 - t2 

gi = , 92 = , 93 ^ - — — = +1 , 

i — 12 

SO the boundary (b) parameterizes one single point; 

(c) h G [0, +1) , t2 = = 
91 = , 92 = , 93 = +9 = *^ 



2-ii ' 

(d) h e [0,+l) , t2 ^ -oo ^ 



91 = ^/tl{l - ti) cos (j) , 
92 = A/ii(l - ii) sin , 93 = 2ti - 1 . 

It is ellipsoid 4(9^ + 92) + (/i = 1 passing through two points 

ii ^ 1 91 ^ , 92 ^- , 93 ^ +1 , 

ti ^ 0, 91 ^ , 92 ^ , 93 ^ -1 • (59) 

The structure of the boundary may be illustrated by Fig. 2. Some clarity can be added with 

the help of the inverse formulas (let qi = q rii): 

9(713 + 1) 9(n3 - 1) , , n2 

h = — T- , t2 = — , tan = — ; 

1 + 9 1 — 9 "1 

from whence it follows 

2q 



n3 = +1 , ii = , t2 



g-1 ' 



29 

m = -1 , ti = — — , i2 = . 
9 + 1 



(d) 





M 


0)-space 




(0) 








(a) 
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-space 




t2 



Fig. 2. Parabolic coordinates in i?3 



In particular, at 7^ ±1 the boundary q = {\ — A), A — t- +0 is parameterized according to 

(l-ns) 



(l + na) 
ti — ^ ) *2 



, tan 



^2 

ni 



2 ' ■ A 

Now we are ready to construct the bound states solutions. In the firs place, note that to have 
solutions vanishing in the origin qi = 0, we must take ai and 02 positive 



ai = +- 



m 



a2 = +- 



m 



(60) 



2 2 

To have solution single- valued and continuous in the region ^ +1 , qi = , 92 = (the 
boundary (b) in the Fig. 2), we must takes positive 61 and negative 62 : 



I -e + e m? 
h^ = +\l—— + — -k^>Q 



-e + -e 



(61) 



besides one should check that the total negative power (02 + 62) compensates a positive power n2 
of the quantum number of the main term of a polynomial at infinity 



a2 + b2 + n2 <0 . 

Additionally we assume positiveness of two expressions under the square roots in (161 p . 
To obtain a polynomial in the variable ti, we require 



(62) 



(63) 



/3i = -m , ni = 0,1,2,3,... 
In turn, the same in variable t2 can be reached by 

h = -n2, 712 = 0,1,2,3,... (64) 
Equations ()63p and (j64p will give (let A'^i = 2ni+ | m | +1 and A'^2 = 2n2+ | m\ +1) 

y2(7^^^yT^^^2T^ = yrT2(+i^ - iVi, 

v/2(-e - e) + m2 - 4A;2 = iVa - .Jl + 2{-e-e), 
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Squaring both ones, after simple manipulation we get 

[ Vl + 2(+e-e) -Nif = [N2- y'l + 2{-e-e) 

so that 

Vl + 2(+e-e) -Ni = N2- y^l + 2{-e-e) 



or 



V^l + 2(e - e) + + 2(-e - e) = 2k , ^ = — ^ = ni + n2+ | m | +1. 
Squaring the above equation, we arrive at a quadratic equation (let x = 1 — 2e) 

x2 - = {2k'^ - xf , 

with solution 

g2 - 1 

e = -^ 2 — ' ^ = ni + n2+ I m I +1 = 1,2,3, ... (65) 

The energy levels belong to the interval 

-y<e<(^-e),A:<V^. (66) 

It is a matter os simple calculations to find expressions for ki,k2 and involved parameters. To 
this end, first let us derive simple expressions for roots 



then 



1 e — e 1,, e, /I — e — e 1,, e, , , 



lie lie 
ai + 6i + 2 = -"1 + 2 + ' 02 + 62 + 2 = -^2 + 2 - • (68) 



It is easily checked (j62|) : 



a2 + 62 + n2 < =^ 
-n2 - ^ + ^ (fc - |)] + n2 = + ^ (A; - ^) < ; (69) 



which holds if (see (j66p ) 



fe--^<0 ^ fc<Ve. (70) 



Now, expressions for ai , a2 are 



C 6 

ai = — ni + {k + -) = n2+ | m | +1 + - , 

6 6 

a2 = — n2 + {k — -r) — Ni = ni+ | m | +1 — -. 



And finally, for bi , 62 



61 = ^2(6 - e) + m2 - 4A:i = ^ [ + ^1 + 2(e - e) - iVi 



1 
2 L 



(A; + -)-iVi 



1 

2 

1 



+(n2 — ni) + 



k 



>0, 



62 = "2 v'2(-e - e) + m2 - 4A;2 = - [ +^1 + 2(-e - e) - iV2 



1 
2 



(fe--)-iV2 



-(ra2 - ni) - - 



< 



it must hold restriction n2 > n^. Besides, note that 



61 + 62 = . 



Additionally, expressions for /3i , P2 can be verified 



m I 1 1 
~ 2 2 L 



{k + -)-N, 



m 



1 1 
2 ''"2^2 



{k 



Pi 
1 



ml 1 
V + 2+^^ 



1 -e + e 
4^ 2 



m I +1 - iVi 



-ni 



iV2 



P2 



ml 1 
V+2+^^ 



1 — e — e 
4 + ^ 



m I +1 - 



-n2 . 



Let us find expressions for separation constants ki, k2 



- 2e + 2e + - 4A;i = [ (n2 - ni) + - , 
-2e - 2e + - 4A;2 = [ (n2 - ^i) + ^ ]^ ; 

from whence it follow 



4A;i = (A: + -)2 - [(n2 - ni) + -f + - 1 , 
Ak2 = {k- If - [{n2 - ni) + lf + m''-l. 



In particular, one obtains ki — k2 = —ie 
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VI. THE HYDROGEN ATOM IN SPACE S3, BOUND STATES IN COMPLEX 

PARABOLIC COORDINATES 

The task consists in separating from all solutions found in Section 4 those which describe 
possible bound states. We will need some details in parametrization of the model H3 



i yi = V-hh cos , iy2 = V-hh sin (/) , 

tl + t2 — 2tit2 2 — ti — t2 

^ ys = „ , yo 



2V(l-ti)(l-t2) ' 2V(l-ti)(l-t2) ' 

ti = {y + 2/3) {y + iyo) = a e*"* , t2 = {y- yi){y - iyo) = b e"*"* , 

tan<P=^ , tit2 = yl+ yl . (74) 

y\ 

Let us consider peculiarities of the parametrization with (ti,t2i0) - {a,b,a;(j)). The origin is 
described by 

ti = , t2 = (a = , 6 = 0) =^ = , yo = +1 • 

The closed line t/q + y| = 1 consists of two parts: 

ti = , t2 / (a = , 6 / 0) =^ iy3 = ^^7=^ , yo = ^^=j^ 5 (^5) 

and 

ti / , t2 = 0, (a / , 6 = 0) =^ iy^ = ^^7=^ , yo = ^^^= 5 (^6) 
here the coordinate (/> is a "dump" one. For another closed line 

2/3 = 0, yo = 0, yl + yl = l; (77) 

we have description 



a + 6 = 2 , ti = e*" ^ 1 , t2 = e"" ^ 1 

— 5" 1 , 2/1 = cos , 2/2 = sin 
ti -\- 12 — 2 i , i , 

^"^= 2V2-ti-t2 = 2^'"^^"^^ = 71^'"'""^° 
2 - ti - t2 1 , 1 



This means that to the line (j77p there corresponds one single point (1, 1) on the boundary (see 
Fig. 1), and here the coordinate (j) is not "dump". There exists one other peculiar region to which 
the boundary a + 6 = 2 (see Fig. 1) is referred 

2/0 = 0, 2/? + 2/i + 2/i = 1 , (79) 
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at this 

sin a = , cos a = 1 , ti = a , t2 = b =^ 
ti + t2 = a + 6 = 2 , </> G [0, 27r] , 
yi = +v^a(2 - a) cos<f> , 
y2 = +\/a(2 - a) sine/) , = a - 1 . (80) 

Finally, we must remember on special additional restriction to which obey two complex coordi- 
nates 

(1 - t2) = (1 - h) (-|) , (1 - h) = (1 - t2) (-|) . (81) 

In particular, these means that if ti — t- 1 it then t2 — ^ 1 T 0, and inversely. 

Now let us separate solutions for bound states. To have functions vanishing on the axis yi = 
0,2/2 = 0, we must take positive ai and 02: 

I "I I I I /^^N 

ai = +^ , as = +^ . (82) 

Let us consider behavior of the functions on the line yf + t/l = 1 ^ see ([78]) . Because here 
ti — )• 1 , ^2 — ^ li to have continuous and finite solutions we must impose restrictions 



Re 61 > , 61 = + ^/ + — -ki 



e — ie m? 



€ + ie m? 



Re 62 > , 62 = + Y + — - fc2 . (83) 

In general, instead it is enough to require only 

Re (61 + 62) > . (84) 

Assuming 



Re \ + >0, Re +Ji + i±^ >0, 



let us reduce hypergeometric functions into polynomials 



/?i = -rii , ni = 0,1,2,...; 

h = -n2, ns = 0,1,2,...; (85) 
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thereby it is supposed that there exist such values for energy and complex ki , k2 at which imaginary 
parts of two square roots in /3i and f32 will cancel out each other. Eqs. (i85]l give 



I ml e — ie m? , 1 /l e — ie 

01 = + \ \ ki-\ a/-H = -ni , 

^ 2V24 2V42 



or differently 



I m I e + ie m? , 1 /l e + ie , . 



m I +1 + 2ni + V2(e - ie) + - Aki = +^Jl + 2(e - ie) 



m 



+l + 2n2 + y2(7TiiyT^^^234fc^ = +yrq:^(7^~) . (87) 



Let 



Ni = 2ni+ I m I +1 , iV2 = 2n2+ | m | +1 , 
then eqs. (f87|) take the form 

+ ^J2e - 2ie + - 4ki = +y^l + 2{e-ie) - Ni , 



+ ^2e + 2ie + m? - 4^2 = + 2(e + ie) - N2 . (88) 
Squaring both equations, after simple manipulation we get (remembering on ki — k2 = —ie) 

[ +^/l + 2{e-ie) -Nif = [ +^/l + 2{e + ie) - A^2 ]^ 
from this two equations follow 

Vl + 2(e - ie) -Ni = I + 2{e + ie) - N2] , 



Vl + 2(e-ie) -Ni = + 2{e + ie) - N2] 

They give respectively 



x/l + 2(e - ie) - Vl + 2(e + ie) = Ni - N2 = 2n 



Vr+2(r^ + Vl + 2(e + «e) = Ni + N2 = 2k . (89) 

Note, that the first of them cannot be valid, further we will consider only second equation in 
89|) - it represent a correct rule for energy quantization. Squaring it we obtain 



(1 + 2e) + V(l + 2e)2 + 4e2 = 2k'^ . (90) 
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Its solution is (let 2e + 1 = x) 



+ 



1 



k = ni + n2+ I m I +1 . 



2P ■ 2 

With the used of the formula for energy levels 

2e + l = + k'^ , 
k-^ 

one can derives simple expressions for all involved parameters. First, one finds 





/l e — ie 




'4+ 2 






/l e + ie 


+\ 




14^ 2 



1 



. e , 



+- \/^^ ~ -p ~ = +- {k — i-) , 



k' 



1 



Now expressions for 



and then ai,a2'- 



1 1 ■e^ 

ai + bi + - = -ni + 2^ k^' 

,1 1 .e, 

02 + 02 + 2 = -»^2 + 2 *• ^ ^^-^ ' 



.e II-, .e 

ai = —ni + [k — I—) = 722+ m +1 — z 



0:2 



For 61 , 62 we easily produce 



1 



k ' 

6 C 

n2 + {k + i—) = n\+ \ m \ +1 + i— . 



k' 



1 



k 



1 



61 = +- V2e - 2ie + rin?-Aki = -[ +^Jl + 2(e -ie)-Ni] = - +{n2-n{)-i- 

A Zi Zi \- iZ 



1 



1 



1 



h2 = +- + 2ie + m2 - 4fe2 = - [ + V^l + 2(e + ie) - N2 ] 
It should be noted that an identity holds 

61 + 62 = ; 



(n2 - ni) + i- 



(91) 



(92) 



the latter is enough to ensure finiteness of the solutions at the region ti ^ 1, f2 ^ 1- 

It was emphasized above that it has sense to examine continuity properties in accordance with 
(see Fig. 1) 

= f2 = a + 6 = 2 =^ (1 - ii) = -(1 - t2) ; 

from whence if follows that if (1 — t\) = 0, then (1 — ^2) = 0. In other words, we have no ground 
to expect continuity of the following types 



*(ti^l,t2), or ^{h,t2^1). 
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Additionally, one can check expressions for /3i , /32 



ml 1 1 
V + 2 + 2L 



1 1 

2 2 



(/c + i- 



m I 1 
V + 2+^^ 



1 e — ie 
4 + ^ 



m I +1 - iVi 



-ni, 



ml 1 
V + 2+^^ 



1 e + ie 
4 ^ 2 



1 , .e I m I +1- ATa 



Finally, it is the matter of simple calculations to specify ki and ^2 : 



2e - 2ie + - 4fei = [ +(n2 - ni) - ]^ , 



2e + 2ie + - 4A;2 = [ -(n2 - ni) + i- ]^ 



that is 



4ki = {k- iyf - [ (n2 - ni) - ]2 + - 1 , 
4/c2 = {k + i^f - [ (ra2 - m) - ]^ + - 1 . 

In particular, the identity ki — k2 = —ie holds. 

VII. THE RUNGE-LENZ VECTOR AND PARABOLIC COORDINATES 

At separating the variables in Schrodinger equation two constants were introduced ki ^2; the 
problem is to find an operator that is diagonalized on wave functions (31) with eigenvalues {ki + k2) 



B h /2 e*'"'^ = (fei + k2) h h e 



im<j> 



(93) 



For the operator^ one can obtain the following representation 



^ /, X ^ ^ 9 (-H + ie) 1 52 

^ d , , d i-H-ie) 1 52 

or after substituting the expression for H 

^ . h+t2-2ht2 ^2t2{l-h){l-2ti) d 2ti{l - t2){l - 2t2 
JD = — le ■_ ■_ 1 ] ] — — h 



d 



ti-t2 



tl-t2 



dti 



^ 2tM^-ti? SP_ ^ 2t2h{l - t2f 52 



t2 - ti dt2 
tl+t2 ^2 



tl-t2 



ml 



t2-tl 



dtl 



2tit2 



(94) 



(95) 
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note the identity 

. ti + t2- 2tit2 . „ . q3 

— le = — ze cosfc' = —le — . 

h-t2 q 

Now we turn to establishing connection between B and Runge - Lenz vector. It is convenient 
to solve the task in the same time both in space and S3. 

In and S3 the quantum mechanical Runge-Lenz operator is constructed from momentum 



and orbital momentum by the formula 



where 



A = el + U[LP]-[PL]), (96) 

Pi = T Qiqj) ^ , L=[qP] , (97) 

upper sign corresponds to the model H3, lower corresponds to 5*3 model; operators L and P are 
measured in units h and h/p respectively. 

In correspondence with symmetry of space models, the components of P, L obey the commu- 
tation relations of Lie algebras so(3.1) and so(4): 

[La, Lb] = i eabc Lc , [La, Pb] = i eabc Pc , [Pa, Pb] = ±i eabc Lc . (98) 

Since in the above expression for B specific term —ie qs/q is present, (in the model H3 we see 
the term e q^/q), it is natural to look for certain relationship between B and ^3. 
Rather long calculation give the following result 

in H3, B = {A + P) , in S3, iB = {A + iP) . (99) 

VIII. DISCUSSION 

In should be emphasized that the possibility to employ complex-valued coordinates in space of 
positive constant curvature 5*3 can be used in other coordinate systems as well - Olevsky's results 
provide us with 34-6 = 28 such special cases. For instance, a complex analogue in spherical space 
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58l |. ant it was used 



5*3 for horospherical coordinates of Lobachevsky space was introduced in 
to examine Shapiro's plane wave solutions of the Schrodinger equation in spaces 5*3 by analogy 
with H3 model. 

Such a possibility evidently will extend the class of integrable problems in these spaces - see 



53] • Also, complex coordinates in 3D-spaces of constant curvature can be of interest in the context 



of the theory of the Lorentz group 50(3, 1) - see in 



5S 
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